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A new diagrammatic representation for correlation functions in the in-in formalism 
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In this paper we provide an alternative method to compute correlation functions in the in-in 
formalism. We develop a modified set of Feynman rules that allows for an easier computation of loop 
corrections. With respect to the standard in-in formalism (where the diagrammatic representation 
is very compact but the evaluation of each diagram is rather involved and contains many different 
contributions), in our formulation the major effort is devoted to drawing the diagrams, while their 
mathematical interpretation is straightforward. Our method is particularly suitable for applications 
to cosmology. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

There has been an increasing interest in recent years 
towards the study of non-linearities during the inflation- 
ary era. It is very well known that inflationary theories 
predict the behavior of the inflaton and metric fluctu- 
ations to be fairly approximated by that of very light 
fields with nearly Gaussian correlation functions. Still, 
although severely constrained by observations [l|, the 
non-Gaussianity of the primordial cosmological pertur- 
bations is not ruled out. Some level of non-Gaussianity 
is indeed expected, as a consequence of the intrinsic non- 
linearity of the theory at least in the gravitational part of 
the action. However, corrections to the linear behavior 
are suppressed by the smallness of the couplings, which 
involve the Newton constant and the slow-roll parame- 
ters Hid, [I]- This makes their predicted value be so small 
that the detection of non-Gaussianity would be very dif- 
ficult, at least for the simplest models of inflation when 
just one scalar field is involved. 

Recently, it has been suggested Hf that larger non- 
linear contributions to the statistics of the cosmological 
perturbations could in principle come from loop correc- 
tions to the correlation functions. The interest of study- 
ing loop contributions is that they contain time integrals 
over the total duration of inflation, and might thus have 
a dependence on the scale factor possibly compensating 
for the smallness of the coupling constants. Actually, 
when contributions beyond the tree level are taken into 
account, more powers of the coupling constants are in- 
cluded with the inclusion of additional interaction ver- 
tices; nevertheless, their increase might be compensated 
by the time integrals in the loops. At the end of in- 
flation, loop contributions could in principle grow large, 
modifying the tree-level behavior of the correlation func- 
tions. Even though such integrals over time are not able 
to produce positive powers of the scale factor aij) (which 
would easily overcome the powers of the coupling con- 
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stants), but can grow at most as powers of loga(r), some 
work is still needed in that direction. 

Unfortunately, the mathematical aspect of computing 
quantum loop corrections in cosmology is more involved 
than in ordinary quantum field theory. The reason is 
that the correlation function we want to compute are 
true expectation values between the same in vacuum and 
not in-out transition amplitude. Actually, the very same 
concept of an out vacuum is misleading in cosmology, 
since there is not a moment when one can safely consider 
the interaction to be switched off. 

When one wants to compute in-in expectation values, 
the Schwinger-Keldysh formalism 0, [1] has to be used. 
In this formalism, the time is thought to evolve along a 
close path, from the remote past to some finite time (af- 
ter the time at which expectation values have to be eval- 
uated) and back. The evaluation of Feynman diagrams 
involves an extended set of rules, since two different fields 
and 4>- are introduced in the forward and backward 
time branches respectively. As a consequence, one must 
distinguish positive and negative vertices describing 0+ 
and (p- interactions. The propagator between two posi- 
tive vertices is time ordered as usual, while the one con- 
necting negative vertices is anti-time ordered. Finally, 
mixed propagators also exist, in which by definition 
is following (pj^. in time. This fairly complicated set of 
rules makes the computation of the expectation values 
rather cumbersome, due to the presence of many more 
terms than for the in-out transition amplitudes in stan- 
dard QFT. 

A different formulation of the in-in formalism was 
given in This approach is particularly suitable for 
applications to cosmology because it automatically takes 
care of several cancellations in the late time behavior of 
the different contributions, which are not immediately 
apparent in the usual formulation. This is achieved writ- 
ing each term of the perturbative series as a product of ei- 
ther commutators of two free fields or expectation values 
of two free fields. Unfortunately, lengthy calculations are 
still needed in order to obtain these reformulated contri- 
butions to the correlation functions. Although straight- 
forward, this procedure can be very long and involved. 
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especially for high order terms. 

The aim of the present paper is to provide a graphi- 
cal representation reconnecting this formulation to some 
modified Feynman diagrams. The most tedious part of 
the calculations can be done graphically, thereby combin- 
ing the useful features of the new formulation with the 
power and efRciency of the diagrammatic expansion. 

The plan of the paper is the following. In Section |TI] 
we review the basics of the in-in formalism. In Section 
mil we develop an alternative formulation leading to the 
same results. In Section IIVI we provide a graphical rep- 
resentation of our procedure, and in Section |V] we make 
the connection with a modified version of Feynman's di- 
agrams, whose symmetry factor is discussed in Section 
IVII Finally, we draw our conclusion in Section I VIII 



II. THE STANDARD FORMALISM 

We consider the action for a self-interacting scalar field 
in a fixed de Sitter space-time background. The metric 
for such space-time is 



As' 



dt^ - a2(t)dx2 



(1) 



where the scale factor a{t) evolves exponentially like c^*, 
so that the Hubble parameter a/a = H is constant. This 
action reads 



dxa^ 



^m'^^ - V{cf,) 



(2) 



giving, when variatcd with respect to (j), the equation of 
motion 



(3) 



The retarded Green function for the Klein Gordon 
equation is the solution of the inhomogeneous equation 

df + 3Hdt-^+mAGn{x,x')^^-^^^; (4) 

it is given by 

Gfl(x, x') = i^{t - t') [(/.o(x), 0o(2:')] , (5) 

where the square brackets represent the commutator 
(j)Q{x)(l)Q{x') — 4'o{x')(j)o{x) and (jiQ is the canonically quan- 
tized free field: 



h[x) -/dk [e''""0fe(i)ak + h., 



(6) 



where the modes 0k (i) are the solution of the linear equa- 
tion (obtained with V = {)) 



4)k + iH(t>k 



m'(f>k 



0. 



(7) 



The objects we want to calculate are the in-in cor- 
relation functions at equal times. As shown in the 



expectation value between two in vacua of a product of 
r fields Q{t) = (f>{xi) . . . 4>{xr) all evaluated at the same 
time t ^ ti = ■■■ = tr is given by 



N=0 •'^'"^ 

X (\Hi{tN),-- 

oo 



dt 



N 



tin -^tin 

Hi{t2),\Hiih),Qo{t) 



(8) 



N=0 



in this equation, Qo{t) = (f)o{xi) . . . (f)o{xr) is the product 
of the same m fields in the interaction picture, which was 
indicated as Qi{t) in and Hj{ti) is the interaction 
Hamiltonian. 

Let us assume for simplicity that the interaction 
Hamiltonian is 



HiiU) 



X 



dya^t,)(b^{U,y); 



(9) 



in this case, the A^-th term of the series ([5]) will be the 
expectation value of 



QNit) 



iX 



N N 



Y[ [dy^a\u)l}{U-l-U)] 



00 (yA^ ) , • ■ ■ [-^0 (y2 ) , [0^' (2/1 ) , Qo (i )] ] ] , (10) 



where we have defined to = t and dyi = dtidyi. In 
particular, choosing Q{t) = 4>{x) we can write 



'^(^) = E hr / n [dy. «'(^.) m-i - u)] 

N=0 1=1 

X ['/>?(2/Ar),---[0o(y2),[0^(yi),0o(.T) 

00 

= 4>i)[x) + E Hn{x). 



(11) 



JV=1 



Since each interaction is proportional to the coupling 
constant, the contribution of A'^-th order in A, i.e. Q^it), 
contains N nested commutators. All these commuta- 
tors involve multiple products of fields, and should be 
transformed into simpler commutators of two fields. The 
reason for doing this is twofold: first, as shown in Q 
commutators of two fields have a different late time be- 
havior than fields themselves, and it is useful to compute 
how many objects of each kind there are in each contribu- 
tion. Second, commutators of two fields arc just numbers 
commuting with everything, and can thus be pulled out 
of any further commutator and expectation value. After 
all the nested commutators are worked out, Eq. pUj) will 
split into the sum of several terms, each of which contains 
the product of N two-field commutators times the expec- 
tation value of the remaining free fields. This calculation 
can rapidly become rather long and involved. Our goal 
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is then to find a simple method to figure out which com- 
mutators and which expectation values of fields should 
appear in each term, with no need to go through all the 
commutator algebra. 

We will now show how the iV-th order term (5jv(0 
of the r-point correlation function can be expressed as 
a combination of the S(j)i^s appearing in the expansion 
([TT|l of the field. Without loss of generality, we can set 
r — 2 and consider the two-point correlation function 
{(j){xi)(j){x2)) ■ The generalization to a higher number of 
fields is straightforward. 

The 1-st order term contains the commutator 

[4>o{y).Mxi)Mx2)] = Mxi) [0o(2/),0o(a;2)] 

+ [ct>l{y),Mxi)]Mx2); (12) 

after being expressed in terms of commutators involving 
only one external field at time and plugged back into 
([Tojl . it yields 



(13) 



Ql(^) = 4'o{xi)5(l)i{x2) + 5(l)i{xi)4>(){x2). 

The second order term involves the commutator 

= Mxim{y2),mvi).Mx2)]] 

+ my,lUx2)m{y2),Ux,)] 
+ my,),Ux,)m[y2),Ux2)] 
+ [(/.Jf(2/2),['/'o (yi), 0o(a;i)]] Mx2); 



(14) 



when inserted into ()10p . the first and last terms immedi- 
ately give 4>o{xi)54)2{x2) and 54)2(xi)4>q(x2) respectively. 
The integration variables yi and y2 in the third term can 
be freely exchanged, which would make the second and 
third term equal. However, the exchange also modifies 
the respective step functions, whose sum now gives 

d{t - ti)d{h - t2) + ^{t - t2)^{t2 - h) 

^d{t~h)d{t-t2); (15) 
therefore, we obtain 



hit) = (i>Q{xi)5(t>2{x2) + <^(?!>2(ii)0o(a;2) 
4- 5(i)i{xi)5(f)i{x2) ■ 



(16) 



This very same procedure can be iterated up to arbi- 
trary order, yielding 



N 



QN{t) = 54>j {xi)6(j)N-j {X2), 
3=0 



(17) 



where the 0-th order terms are meant to be S(j)oixi) = 
4'o{xi) = ip and S(f)Q{x2) = 00 (2^2) = X- Summing up all 
perturbative orders one gets 



Q{t) - }^ QNit) 

N=0 



^(5(/),(.Tl) 



2=0 



^S(t>j{x2) 



(18) 



which is what we could expect since we are considering 
Q{t)=(l){xi)cj){x2). 

Generalizing to the case of r fields one has 

(Xi) ■ ■ ■ d(t)j^{Xr), (19) 

\j\=N 

where the sum is taken over all the vectors j £ N'' such 
that \j\ = ji + ■ ■ ■ + jr ~ N. Again, this gives 



Q{t) = Q^(^) = n 

7V=0 j=0 



1=0 



(20) 



It is thus clear that, in order to reconstruct the expan- 
sion in powers of the coupling constant A of any r-point 
correlation function {(j){xi) ■ ■ ■ (j){xr)), one can just com- 
pute the expansion (|11[) . Then one simply has to take 
the product of all the possible combinations of the (50i's 
having the same overall number of powers of A. There- 
fore, our goal will be that of giving a simplified way to 
calculate (fTTj) . 



III. PERTURBATIVE EXPANSION OF THE 
EQUATION OF MOTION 

We now want to find a simple method to calculate 
the evolution of the interacting field as an expansion in 
powers of the coupling constant, in order to have a for- 
mulation equivalent to (jlip but still not relying on the 
involved computation of all the nested commutators. In 
this Section we show that this can be achieved simply 
with a perturbative solution of the equation of motion, 
once the order of the non-commuting pieces of the ex- 
pansion has been appropriately taken care of. 

Let us then look at the same problem from a differ- 
ent perspective, and solve perturbatively the equation of 
motion ([3]) for the operator field 



4){x) = 4)o{x) + S(j)i{x) + S4i2{x) 



(21) 



expanding around the solution with A = (the field in 
the interaction picture). Higher order contributions can 
be obtained by iteration, plugging the known solution 
back into The next order contribution will then be 
the solution of a linear differential equation with a source. 

The evolution equation for S(f>i (considering for sim- 
plicity a cubic potential V = ^0'^) is 



A 

'2' 



(22) 



whose solution can be found with Green's method. Using 
the retarded Green function (l5|), one gets 



(x) = I dyat,GB{x,y)^t{v) 



(23) 



When one wants to find the second order correction 5(j)2 , 
the potential term in ^ becomes V'{4)o + (50i) and must 
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be expanded up to second order in A. In this case one 
must pay attention to the fact that (po and S(pi do not 
commute, since they are quantum operators and contain 
free fields evahiated at different times. The expansion 
yields 

V'{cj)o + 5^1) = [-^0 + <^o¥i + + 0(A')] ; (24) 

keeping only the second order terms one has the source 
for the equation of motion for 5(j)2^ which reads 



m 5(f)2 



X 



{0o,<^<^i}, (25) 



where the anticommutator |(/)o, = (j>oS(pi + Scfii (j)o 
symmetrizes the dependence on the non-commuting 
fields. The solution can be found again with Green's 
method, and using the expression (|23p for Scjji gives 



6(j)2{x) 



dya^ I dzal 



X GRix,y)GRiy,z){My).^liz)}, (26) 

Similarly, from the third order contribution to the poten- 
tial one gets 



J dyalGR{x,y)[{My),SMy)}+Sciyl{y) 
[Mv)^ Gniy, ^){0o(2:), Gi^(z, w)4>l{w)}j 



+ GR{y, z)(j)l{z)GR{y, w)^o(w') 



(27) 



More generally, including corrections to the field value 
up to A^-th order, the potential term in ([3]) can be written 
as 



N 



A 



N 



,4=0 



N 

E 

j=i+i 



, (28) 



which allows to express iteratively the general solution 
as 



-^Gr 



>-E 

i=o 



E 

J=0 



j,5(t>2i-j} 
^,"^21+1-^} • 



(29a) 
(29b) 



In general, any contribution to (jlip is an integral of 
a retarded Green function followed by two lower order 
S(f>i's. If these are not both (j)Q, the first Gr is followed 
by another Gr, and so on. The generic contribution to 



6(I)N thus contains a tree of A'' nested Green functions 
stemming from the first one. Each of them is followed 
by a pair of objects, which can be free fields or other 
Green functions. It also contains N integrations over the 
N vertices of the tree. When the two objects after a 
Green function are not symmetric (like a free field and a 
Green function, or two Green functions followed by dif- 
ferent objects), they must be symmetrized with an anti- 
commutator. The total 64>isf is given by the sum of such 
contributions over all possible inequivalent dispositions 
T/v of the N Green functions. 

When a generic interaction (/)" is considered, the pre- 
vious results still hold. The only difference is that any 
Green function of the tree is now followed by n — 1 ob- 
jects. The symmetrization of asymmetric objects must 
now be done by mean of a generalized anticommutator. 
For n — 2 (out of n — I) equal objects, this is defined as 



in— 2 in—1 
01 02 



51(^2 



in — 3 in— 1 1 
h 02 01 



■i.n— 1 171 — 2 
02 01 ' 



(30) 



which reduces to the standard anticommutator for n = 3. 
In an arbitrary situation with ki free fields, k2 groups of 
n — 1 free fields and in general ki groups of i{n — 2) -I- I 
free fields (with fci -I- ^2 + • • • = ^t- — 1) , this generalized an- 
ticommutator gives the (n — l)!/(fci!fc2! • • ■ ) inequivalent 
configurations. The generic (50jv then looks like 



S(f>N{x) 



rr ^ 

^) E/ Il'^y^-'ii^) 

/ _ L A — 1 



GRix,yi) 

x{...,GRiyN-i,yNW^-\yN)}, (31) 



where a number from to — 2 of (generalized) anticom- 
mutators is needed depending on the number of vertices 
that has to be symmetrized. The total number of free 
fields in each term of the sum is A^(n — 2) -I- 1, since each 
vertex of the tree one of them is replaced by a Green 
function carrying n — 1 more free fields. 

One can check that the iterative solution of the equa- 
tion of motion gives the same results as the expansion 
PT|) . after all the nested commutators have been re- 
duced to products of commutators of two fields. In or- 
der to shorten the notation, we set 0o(a^) = ip, and 
4>o{yi) = ipi (not to be confused with the i50i's, for 
which the index stands for the order of the expansion!). 
The first order term in (|lf[) contains the commutator 
[ipi , -0] = 3-01 [ipi , ■0] ; recalling the definition ([5]) of the 
Green function, the complete integrand reads 

^i9(t-fi) [03(2/1), 0o(a:)] =-^Gfl(x,yi)02(yi); (32) 



once integrated, this is identical to the result 
obtained with the expansion of the equation of mo- 
tion. Similarly, the second order term of (|f ip contains 
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[V-i, = 32{V'i,V'i}[V'2,i/'i][V'i,V'], and yields 



A 



^[^ 2) 'fi) {'^0(2/1)''^^ (2/1 ' 2^2 )'/>o(y2)}; (33) 

here again each two-field commutator times the appropri- 
ate step function has been turned into a Green function, 
returning the same expression as the integrand in (|26p . 

If we want to go to the next order we need to add 
one more nested commutator to those we just com- 
puted. The two-field commutators [^/'2, V'lJlV'ii ^] in 
the inner part are just numbers, and can be pulled 
out of the outer commutator. The latter then acts 
only on Since it obeys the Leibniz rule 

it will act separately on the two sides of the anti- 
commutator, giving {[tpi,ipi],'ipi} + {V'l, [V'i'V'i] } = 

3{ {^1 , 1^2 } , V'l } [V'a , V'2] + 3{ V'i , } [V-s , V'l] • Adding the 
step functions, the constants and the remaining two-field 
commutators, the first term becomes 



A 



Gr{x, yi)Gi?(2/i, ?/2)Gfl(y2, ya) 



x{Myi).{My2),^liy3)}}, (34) 

equalling the first term in (j27|) . In the second term, the 
■iJ's do not immediately match the two-field commutators. 
However, 1/2 and are integration variables and can be 
exchanged, which gives i?(ii — t2)i?(t2 — *3){V'3iV'2} = 
'&{ti — t2)^{ti — is)^'!^'!; the step functions both 
start at ti, and we obtain 

aV 

) GbXx, yi)Gi?,(yi, ^2) 0o(y2)Gfl,(yi, ys) (pHys) 

(35) 

as in the second term of (|27p . 

At each step, a higher order term can be obtained from 
a lower order one by turning one of its free fields into a 
new Green function. As shown, this comes from the fact 
that when the new outer commutator acts on the free 
field it gives a two-field commutator times two new free 
fields. One can show that the temporal step function can 
always be rearranged so that the coordinates of one of 
them match those of the new two-field commutator in 
order to produce a Green function. Finally, the Leibniz 
rule provide the necessary symmetrization. All of the 
above can be applied as well to the case of a generic </>" 
interaction. We therefore conclude that the pcrturbative 
expansion of the equation of motion for the operator field 
reproduces the very same results as the in-in formalism 
in standard QFT. 



IV. GRAPHICAL REPRESENTATION 

As one can see, solving the equation of motion at 
some given order rapidly produces a very high number 



of terms, and also this procedure can thus become very 
complicated. Therefore, we would like to find an auto- 
mated method to reproduce efficiently all contributions 
to each term of (jlip , and reconstruct their mathematical 
content. In this Section, we introduce a graphical rep- 
resentation describing the pcrturbative expansion pip . 
This will reach our goal, allowing to construct pcrturba- 
tive corrections to the field evolution up to any order in 
A simply as a combination of graphs. 

The field in the interaction picture, and the first and 
second order correction ((23|) and ([26]) can be represented 
respectively as 



(l)o{x) 
S(l>i{x) 



X y, 



(36) 
(37) 



and 



5^2(2;) 




(38) 



where the rules to reassign to each element of the graph 
its original mathematical meaning arc: 



1. a dot is associated with the field coordinate 



x; 



2. each solid line represents a retarded Green function, 
the farther end from the dot being in the past: 



X Vi 



Gfl(x,2/0, ^—^GR{y,,y,)- 



3. each vertex carries an integration factor: 

4. each dashed line with a cross carries a free field: 

---X = 0o(2/i) ■ 

For example, the interpretation of the graph p7p giv- 
ing the first order correction is 



A 



GR{x,y) 



(39) 



'/'o(y) 



in this case, the two free fields are evaluated at the same 
time, and their overall contribution is simply 4'oiy)- This 
correctly reproduce the result obtained in 
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As for the second order contributions, the first of the 
two graphs in ([38|) gives 
, 2 

^3 J — 3 




0o(y) 



(40) 



where we get a double integration factor due to the pres- 
ence of two vertices, and two retarded Green functions. 
In this second order graph one also has to care about 
how to order the free fields, since these are now evalu- 
ated at unequal times and do not commute: we could get 
either </'o(2/2)0o(j/i) or (/'o(2/i)0o(j/2)- Choosing to order 
the fields "from top to bottom" would give 




(41) 



the second graph in p8|) can be obtained from the pre- 
vious one simply by flipping the lines in the first bifur- 
cation, and thus differs only for the order of the fields. 
With the same convention we would have 




(42) 



and the sum of the two graphs correctly reproduces d(j)2 
as given in p6)) . 

Choosing a different convention to order the fields (for 
example, "from bottom to top" ) would reverse the contri- 
bution of each graph, but their sum would be unchanged. 
Therefore, independently of the rule chosen for the single 
graph, the sum of the two mirroring graphs will always 
give the anticommutator. We also notice that the sec- 
ond order graph contains two bifurcations, but only one 
of them (the one labeled by y in our graphs) introduces 
an anticommutator. The second bifurcation is symmetric 
(like the one in the first order graph) and does not have 
any effect on the field order when fiipped. In general, 
only bifurcations with distinguishable branches change 
the graph when they are flipped and introduce anticom- 
mutators. It would be useful to denote the sum of such 
graphs as 



X y(^^-Y. 

\ 



(43) 



where the index stands for the number of graphs pro- 
duced by the anticommuting bifurcation. Independently 



of the convention chosen to order the free fields, this 
symmetrized graph directly yields the anticommutator 

{My)Al{z)\. 

In general, each correction S(j)N(x) is represented by a 
tree-shaped graph with N bifurcations connected by lines 
representing the retarded Green functions. Higher bifur- 
cations therefore correspond to points in the past light 
cone of the space-time coordinate x and of lower bifur- 
cations. The final branches of each tree graph represents 
all the free fields contained in the expansion, of which 
we will take the expectation value. The previous exam- 
ples suggest that asymmetric graphs should always be 
considered together with its mirror counterpart (s): the 
sum of all the graphs lead to the correct symmetriza- 
tion procedure of these non-commuting free fields. More 
precisely, one has to consider the sum of all the graphs 
that can be obtained from each other simply by flipping 
the branches of one or more bifurcations, the sum of 
all graphs will then contain as many anticommutators as 
distinguishable bifurcations, with each anticommutator 
involving the fields (or groups of fields) attached to the 
branches that are being flipped. 

Following the same steps we can write the third order 
expansion (|27p as 



Sh{x) 




(44) 



where the double index in the second graph means that 
we now have two bifurcations {y and z) with two pos- 
sibilities each, for a total of four graphs that we must 
sum. Flipping the bifurcation in z exchanges 4>o{z) 
with (^q(u'), while flipping the first one exchanges 4'oiy) 
with all the other fields taken together. This will 
therefore results in the double nested anticommutator 
{(/)o(2/), {0o(z), 0o(^)}}' prescribed by ([27]) . On the 
other side the first graph does not contain any distin- 
guishable bifurcation: indeed, its free field content is 
globally symmetric because, although (j)l{z) and (I)q{w) 
do not commute, z and w are dummy integration vari- 
ables and can be safely exchanged without affecting the 
overall result. 

When dealing with a generic interaction such as </>", 
each vertex will contain n — 1 outgoing branches. If 
only one branch has one more bifurcation attached, then 
there will be n — 1 possibilities of twisting the "(n — 1)- 
furcation" . If the branches containing further pieces 
are two, we will have {n — l){n — 2) possibilities if the 
sub-graphs of the two branches are distinguishable and 
(n — l){n — 2)/2 if they are not. In general, the number 
of possible dispositions of k equal objects in n — 1 slots 
is given by the binomial coefficient ("^^). In this case, 
we need to use the generalized anticommutator ([50)1 . 

We can thus complete our set of prescriptions with 
the following rule, which describes the symmetrization 
procedure of the free fields in asymmetric graphs by mean 
of anticommutators: 
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5. the sum of all the mirroring graphs with I distin- 
guishable bifurcations contains I (generalized) an- 
ticommutators exchanging the group of fields at- 
tached to each branch of the bifurcation. 

We can easily show that this graphical method is cor- 
rect up to any order. First of all, we notice that the 
second order graphs in ([38]) can be obtained combining 
together two first order graphs in the two possible posi- 
tions. Also the third order graphs in (|44)) arc obtained 
attaching two first order or one second order graph re- 
spectively on the branches of a first order one. This re- 
flects the structure of the analytic expressions and 
(f26ll for S(j)i and 6(1)2, where there is an overall Green 
function G^. According to our graphic rules, this is ac- 
tually represented by the "handle" of the first bifurca- 
tion, to which one has to attach other lower order pieces. 
Generalizing, since as shown in p9|) all corrections have 
an overall Green function Gn, then the graphical proce- 
dure of obtaining higher order corrections from combin- 
ing lower order terms on the top of a first order graph 
holds at any order. 

However, constructing A'^-th order graphs is actually 
much simpler, since one can simply write down all the 
possible graphs with N bifurcations in all possible posi- 
tions, without caring for lower order graphs. The equiv- 
alence of the two procedures is implicit in the recursive 
structure of Eq. ((^^ for a generic S(pN- The proof is very 
simple: Eq. (|^^ shows that the A^-th order is made up 
by all the possible combinations of lower order graphs 
containing a total amount of bifurcations equal to A^ — f , 
attached on the top of a first order graph. The lower or- 
der sub-graphs arc disposed in all possible positions, as 
a consequence of the presence of the anticommutators. 
Therefore, if the two procedure are equivalent up to the 
(A^ — f )-th order, then they are also up the A^-th. 

So far, we have been describing the graphical method 
to write down in-in expectation values in the case of a 
cubic interaction, where one only has two-ways bifurca- 
tions. However, the procedure is completely general, and 
can be used to describe a 0" theory where ramifications 
with n — 1 branches exist. The steps one needs to follow 
are strictly the same; the only difference is that there are 
now several different ways to flip the branches of each 
(n — l)-furcation, which is a graphical representation of 
the generalized anticommutator pop . For example, the 
second order contribution in a A^^ theory is 



S(f)2{x) 




(45) 



where the sum over the three possible graphs obtained 
from the permutation of the branches produces the gen- 
eralized anticommutator 

0g(2)0g(y) + My)4{^^)My) + ^Uy)4>Uz). (46) 



Summarizing, the rules for constructing an A^-th order 
graph with a generic interaction are 

• start the perturbative tree from the external point, 
connecting any additional vertex to lower vertices 
or to the external point with a solid line 

• each vertex has one ingoing solid line (connecting 
to the lower part of the tree), and n — 1 (solid or 
dashed) outgoing lines 

• outgoing solid lines lead to higher vertices, outgoing 
dashed lines are final. When connecting an addi- 
tional vertex, a dashed line must be changed into a 
solid one. 

• the total number of vertices must equal the order 
of the graph 

One must then take the sum over all the possible dispo- 
sition of the A vertices, and apply the previous rules f 
to 5 in order to interpret the graphs. 



V. CONNECTION WITH FEYNMAN 
DIAGRAMS 

In this Section, we combine several graphical pertur- 
bative solutions in order to recover closed Feynman-like 
diagrams for the correlation functions. Closed graphs 
will therefore represent the expectation value of the free 
field content described by some combination of open tree- 
graphs. 

As we saw in ((T7)) and the product Q{t) of some 
interacting fields can be expanded into the sum of terms 
involving the product of several S(f>i^s. For example, the 
two-point function can be written as 

oo N 

{^{X,)^{X2)) = E E {S<l^ji^i)d<I^N-j{x2)) , (47) 
Ar=0 j=0 

where each expectation value acts on the product of all 
the Gaussian free fields contained in the Scpi's. Using 
Wick's theorem, all these expectation values can be in 
turn decomposed into several products of two-field expec- 
tation values, containing as usual all the possible pairs of 
the free fields. 

We now want to implement this fact into the graph- 
ical representation. From a graphical point of view, an 
m-point correlation function can thus be reproduced by 
putting together m field expansion diagrams (each of 
which is associated to an external point) and connecting 
all their free fields in all possible ways. We will graphi- 
cally describe the contraction of two free fields as 



y 



y 



(48) 



which as we will see can represent both the expectation 
values {(j)Q{y)(j>Q{z)) or {(j)Q(z)(j)Q{y)) . Let us study some 
examples. 
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The simplest case to study is the mean value of the 
field. We immediately see that all the even terms of the 
expansion vanish, since they contain an odd number of 
free fields. We thus have {6(I)n{x)) = for N even. For 
the first order term, we have instead 



(49) 



where the dashed line represents the only available con- 
traction of the two free fields, which is while the 
meaning of all the other graphic elements (solid lines, 
vertices, points) remains unchanged. Since there is only 
one possibility to contract the free fields, the overall fac- 
tor of i coming from the vertex remains. The result is 
therefore 



dya^GRix,y)(l)liy). (50) 



Another useful example is the three-point function 
{(j){xi)(j){x2)(j){x3)) . At lowest (that is, first) order in A, 
this is given by the contraction of one first order and two 
zeroth order terms: 



{4'ixi)4'ix2)(l){x3)) ~ {5(j)i{xi)(j)o{x2)(j)oix3)) 
+ {<l}o{xi)S(j)i{x2)(l)o{x3)) 

+ ('/>o(a;i)(/>o(a;2)<5(/>i(a;3)) , 



(51) 



while all the even terms of the expansion vanish. 

The first term of ([STjl contains the Green func- 
tion Gfl(xi,y) and a total free field content given by 
4'oiy)'t'oix2)4'o{x3) , whose expectation value is 

{(l>o{y)Mx2)(l>o{x3)) = '^{<I^Q{y)Mx2)){<Po{y)(j)o{x3)) 

+ {4>l{y)){M^2)M^^)). (52) 

The last part of this expression would merely give rise to 
the disconnected contribution {54>i{xi)){4>i:i{x2)4>o{x3)), 
which we will neglect. The connected part of the first 
term of ([51]) can thus be represented as 



xiy,' 



X2 



X3 



f X2 



(53) 



a;3 



where the solid line represents as before the retarded 
Green function Gji{xi,y) (y being the coordinate of the 
vertex) , while the dashed lines stand for the contractions 
{</>o(2/)0o(a;2)) and ((/)o(y)0o(a;3)}- Unlike in (gll), the free 
fields arc now evaluated at unequal times and their order 
is relevant. In particular, the free field contribution 4>^(y) 
contained in 54>i(xi) will always appear before 00(2:2) or 
4>q(x3) in any expectation value. 

This has to be repeated for the other two possible con- 
figurations in ()5ip . The complete representation of the 
first order contribution to the three-point function (|5ip 
therefore becomes 



X2l> 



Xl 



-( + 



Xl 



X3^ 



X2f 



X3^ 



X\ 



X2f 



X3* 



(54) 



where in each graph the solid line is attached to a 
different external point. Having said about the first 
graph, in the remaining two the solid lines represent 
the Green functions GR{x2,y) and Gii{x3,y) respec- 
tively. As for the dashed lines, the free field con- 
tent of the last two graphs is (l)o{xi)(pQ{y)(j)o{x3) and 
(j>o{xi)(l)o{x2)4'o{y)- therefore, the dashed lines stand for 
{(j>oixi)(f>oiy)) and (0o(y)<^o(a;3)) in the second graph, 
and finally for (0o(a;i)0o(2/)) and {cl)o{x2)<po{y)) in the 
third. 

From this example we can extrapolate the rule for the 
interpretation of such "external" dashed lines {i.e. con- 
necting different perturbative expansions evaluated at 
different spatial points): the order of the fields in "exter- 
nal" contractions is the same as the respective position 
in the correlation function of the perturbative expansions 
they belong to. Following this rule, the diagrammatic 
representation (|54p of the three-point function correctly 
gives 

{(I){xi)(l){x2)(f>{x3)) c^xjdya^ 

X Gr{xi , y) (00 (2/)'/'o {X2)) (00 (y)0o (-^a)) 

+ {Mxi)4>o{y))GR{x2,y){<t>o{y)M^3)) 

+ {(f>Q{xi)(j)o{y)){(j)a{x2)(t)oiy))GR{x3,y) (55) 

up to terms of order 0{X^^). 

It is also instructive to examine the two-point correla- 
tion function. For example, combining the two first order 
contributions 6(pi{xi) and (50i(a;2) (which give a second 
order contribution) we have 



Xl 



jy," \ z X2 



Xl 



X2 



y- 



(56) 



where the solid lines give the two Green functions 
Gii{xi,y) and Gr{x2, z), and the dashed lines represent 
the double expectation value (0o(2/)0o(-2))^. Also in this 
case, the dashed lines are of the "external" kind, connect- 
ing free fields in different perturbative expansions. The 
order of the fields in the expectation value is thus fixed 
by their position in the correlation function (0(a;i)0(a;2)), 
since 0o(y) belongs to the expansion of 0(a;i) and (j)o(z) 
to that of 0(^2). Finally, each vertex brings a factor of 
^, and since there are two equivalent ways to pair the 
fields this yields an overall symmetry factor of ^ . 

Another second order contribution to the two-point 
function is obtained combining 502 (xi) with 0o(a;2). 
From the two graphs in Eq. ((38|) one gets 



^1 yf 



X2' 



Xl V' 



X2' 



(57) 



where in both these graphs the solid line contribution 
reads Gii{xi,y)Gii{y, z). The line attached to X2 is again 
an external dashed line: since in both cases all the free 
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fields in 502 (a^i) precedes (j)o{x2), this represents in both 
graphs^ the expectation value {(j)o{z)(j)o{x2)) ■ On the 
other hand, the two free fields (^o(y) and (j^oiz) both be- 
long to (502 (3^1 ) and when we pair them we must pay 
attention to their order. If we use the convention "from 
top to bottom" previously applied in (PT|) and then 
this dashed line reads ((f)Q(z)(f)o{y)) in the first graph and 
(0o(y)0o(^)) in the second one. As for the diagrammatic 
representation of S(j)2, using a different convention would 
reverse the order of the fields in the internal dashed line 
of each diagram; still, their sum would be unchanged. 

It will thus prove useful to use a single diagram to de- 
note the sum of the two contributions in (j57p . Deforming 
the lines one gets 



(58) 



(the two graphs are completely equivalent, since they 
both represent the sum of the two terms). We now have 
here two different kind of dashed lines: again an exter- 
nal line (the one attached to X2) connecting two different 
pcrturbative expansions, but also a dashed line between 
points that arc also connected by solid lines (i.e. be- 
longing to the same pcrturbative expansion) . In the lan- 
guage of closed graphs, this can be dubbed an "internal" 
dashed line. The two dashed lines represent the expecta- 
tion value 

{{My), cj)l{z)}M^2)) = 2 {Mz)Mx2)) 

I II I 

X [{My)M^)) + {M^)My))], (59) 



function therefore reads 



(0(a;i)0(x2)) = (0o(a;i)0o(a;2)) 



A2 



dyal dzal 



GR{xi,y)GR{y, z) ({0o(y), ^0(2:)}) ((/>o(^)0o(a;2)) 
+ GR{x2,y)GB,{y,z) {(j)o{xi)(t)o{z)){{(j)a{y),(l)o{z)}) 



+ GB,{xi,y)GR{x2,z){<po{y)(f)a{z)y 



(61) 



up to terms of order 0{X^). 

From the previous analysis we might be tempted to 
say that we can associate to every internal dashed line 
the expectation value of an anticommutator. Unfortu- 
nately this is not always the case, but more complicated 
symmetrization procedures can occur. As an example, 
consider the fourth order graph 



(62) 



contributing to the two-point correlation function. This 
diagram comes from the expectation value 





X2 



(63) 



that is of a subset of d(f>i{xi) containing four terms and 
of 00(2:2). The total free field content in the expectation 
value is 



where the free fields paired through the external dashed 
line arc "pulled out" of the anticommutator, while those 
paired with the internal dashed line are symmetrized. 

Following the same procedure we can calculate the 
last remaining second order term, which combines (l)o{xi) 
with the two graphs for (502 (2^2 )• This gives 




(where again twisting the whole diagram upside-down 
would return an equivalent one). 

In all the previous graphs there are two vertices and 
two ways of placing the dashed lines. This yields then a 
global symmetry factor of i. The two-point correlation 



We are not considering here the possible pairing of <f>o{x2) with 
4>o{y), since this would not give a IPI diagram 



({02(u), IMw), 4>U^)}}Mx2)) , (64) 



where the internal anticommutator is associated with the 
vertex w, the external one with the vertex y, and the 
dashed lines indicate the pairings that we are actually 
considering in ([62]) . 

The peculiarity here is the fact that the graph ([62|l does 
not pair 0o(u') with 4'o{z) (such pair is contained in a dif- 
ferent graph). The only pairings are 0o(u) with 0o(w) or 
with 00 (-z), and are not affected by the order of the latter 
fields. Therefore, the internal anticommutator does not 
contain any dashed line and merely introduces a factor of 
2. The only anticommutator affecting the dashed lines is 
the external one, that reverses both lines simultaneously: 
indeed, the result is 2[(0.u,0„) (020„) -I- {4>u4>w){4>u4>z)]- 
This differs from our naive expectation of the product 
of two anticommutators by the presence of a factor of 2 
and the absence of mixed terms. We also have a factor of 
(1/2)^ from the vertices and a combinatorial factor of 4 
counting the equivalent dispositions of the dashed lines. 
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The net contribution of (l62l) then is 



A J dyaydwa^dza^dua^GR{xi,y)GR{y,w) 

X GB,{w,z)GB{y,u)^ {(l)o{w)(l)o{u)) {(po{z)(po{u)) 
+ {(l)o{u)(l>o{w)) {<po{u)(f)o{z)) {(j)o{z)(i)o{x2)) , (65) 



and as said above it contains only one symmetrization 
involving two internal dashed lines at the same time. 

In general, all intermediate vertices of the trees {i.e. 
vertices with at least one outgoing solid line besides the 
ingoing one) are associated to anticommutators that ex- 
changes the n — 1 slots of the outgoing hnes^. Any dashed 
line connected with both ends to two of the n—1 slots of a 
given vertex z will thus be symmetrized by the anticom- 
mutator. When acting on a solid line, the anticommu- 
tator acts on the whole sub-tree starting from that line: 
therefore not only the dashed lines directly attached to 
the vertex z, but also all those connecting vertices in the 
sub-tree of two outgoing solid lines, will be symmetrized 
at the same time. On the other hand, dashed lines hav- 
ing the second end attached to a lower vertex y (or to 
any other of its outgoing solid lines) are reversed by the 
anticommutator in y and not by the one in z. 

We can now classify the behavior of a generic dashed 
line on the basis of its position in the graph. All closed 
graphs are constructed with open trees of solid lines start- 
ing from each external point, connected with dashed lines 
that close the diagram. As we have seen, the latter can 
be internal or external. Internal dashed lines connect ver- 
tices in the same tree, and are therefore necessarily part 
of a loop with only one dashed line (with all the solid lines 
of the loop belonging to same perturbative tree). We will 
refer to this unique loop as to the symmetrization loop of 
the line. Dashed lines for which such a loop does not ex- 
ist are of the external type, and always connect vertices 
in different perturbative trees. The extreme case of ex- 
ternal dashed line is when the line is directly attached to 
an external point. Any dashed line connecting two points 
X and y is related to the two-field expectation value of 
4>o(x) and (l>oiy), but the order of the fields depends on 
which kind of line we are considering. 

From what said above, it follows that internal dashed 
lines are symmetrized by the anticommutator associated 
to the lowest vertex of their symmetrization loop (i.e. the 
one closest to the external point). Whenever two or more 
dashed lines are symmetrized by the same vertex, the end 
points of all the lines must be exchanged simultaneously. 



^ In open perturbative diagrams not all intermediate vertices are 
symmetrized: vertices with a symmetric sub-tree {i.e. indistin- 
guishable outgoing solid lines) do not bear any anticommutator. 
However, in closed graphs one must sum all the possible ways 
to connect external dashed lines to them, and this provide an 
effective symmetrization of such vertices (see App.|X]for details) 



For example, in the case of ((62|) . the two loops are y-w- 
u-y and y-w-z-u-y and in both cases the lowest vertex 
is y. In general, when k lines are symmetrized by the 
same vertex one needs to introduce the factor (for a (f)^ 
potential) 



i=l 



(66) 



where for convenience we introduced the normalization 
i which will be useful in the following. 

In the case of a generic 0" interaction one has (n—1)! 
permutations of the n — 1 outgoing lines. Any vertex yi 
attached to the line ly- will be associated by each permu- 
tation p e Sn-i to the line p{lyi). The order of two fields 
4>o{yi) ^-iid 00 (-2^1 ) is then reversed when the permutation 
changes the sign oi p{ly-) — p{lz^). The generalization to 
the previous formula to a generic interaction is therefore 
immediate, and reads 

1 ^ ( 



p i—1 



^[p{lz,)^p{ly,)]{Mz^)My^))], (67) 



(n- 1)! 



where again we introduced a normalization factor by 
hand. 

Dashed lines symmetrized by different vertices enter as 
a separate factor, so that there will be as many factors of 
l/(n — 1)! as symmetrizing vertices. On the other hand, 
external dashed lines do not need to be symmetrized: the 
order of the free fields in the expectation value is fixed, 
and is the same as the order the two external points to 
which they are connected have in the complete correla- 
tion function. 



VI. SYMMETRY FACTORS 

We now want to understand how to compute the over- 
all symmetry factor for each closed graph, without going 
through the complicated procedure of connecting all the 
equivalent open diagrams. 

We notice that we assigned in ([66|) the same normaliza- 
tion i to all groups of internal dashed lines symmetrized 
by a single vertex. As explained in App. El this factor 
needs to be introduced when the dashed line connects the 
ends of two symmetric solid lines stemming from a vertex. 
In this case the vertex does not have a proper anticom- 
mutator and the symmetrization of the line(s) is provided 
by the different ways of connecting other lines to its end 
points. The combinatorial factor of the additional lines 
doubles after the symmetrization (the two end points are 
now equivalent), and the ^ is needed to reproduce the 
correct result. When instead the symmetrization is pro- 
duced by a real anticommutator (like in all the examples 
throughout the current Section), the i factor docs not 
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FIG. 1; Fourth order contributions to the two-point correlation function. The symmetry factor is 1 for aU graphs but (5), 
(11), (17), (18), (23) and (24), for which it is |. In graphs (4) and (10) the internal dashed lines share the same symmetrizing 
vertex and are symmetrized together, yielding i [(<^o(m)<^o(w)) {(f>o{u)(f>o{z)) + (</>o(w)(^o('u)) (<^o(z)'?io(ii))], as well as in (17) 
and (23) where they give ^ [((^o(2)0o(«'))^ + ((^o(w)0o(2))^] . In all the other graphs all internal dashed lines are symmetrized 
individually, accounting i [{4>o{yi)4>o{yj)) + {4>o{yj)4>o{yi))] for each line 



arise and must be introduced by hand: the associated 
vertex must then be given a factor of 2 to compensate. 
Finally, vertices with "idle" anticommutators also intro- 
duce a factor of 2. 

It follows that any vertex with just one outgoing solid 
line carries a factor of 2, while vertices with two do 
not carry any additional factor since they are symmetric 
(they could still have of course a combinatorial factor if 
the outgoing lines can be equivalently attached to several 
other vertices, but no factor for the exchange of the lines). 
Final vertices (with only outgoing dashed lines) carry 
the standard combinatorial factor counting the equiva- 
lent dispositions of the dashed lines. All these factors are 
computed holding the solid lines fixed, and must multiply 
the factors of ^ introduced by the vertices, which amount 
to an overall factor of 2~^ for all A^-th order graphs. 

This procedure is equivalent to the one used in stan- 
dard QFT, where a factor of (l/iV!)(l/3!)^ is assigned 
to any graph, and afterward multiplied by the number of 
equivalent configurations. Indeed, the number of possi- 
bilities just for placing the N solid lines (that do not form 
any loop) automatically contains a factor of A^!3^, which 
gives back the original 2^^. The only difference from the 
standard computation of symmetry factors is that one 
has to distinguish between the two types of lines: when 
placing a solid line in a vertex containing also a dashed 
line (and vice versa) one must count the legs of both the 
starting and arriving vertex. This introduces extra com- 
binatorial factors, which exactly account for the factors 
of two we need in asymmetric vertices. 

Let us try to clarify the outlined procedure with a cou- 
ple of examples. In the graph in Eq. (p^ [the same as 
graph (4) of Fig. [1], which starts from (l/4!)(l/3!)'', we 
have 4 • 3 ways to place the xiy line, 3 • 3 for yw (we do 
not count the starting legs in y because yw and yu are 
symmetric) 2 ■ 3 for yu, 3 for X2Z, 2 • 4 for zw (2 starting 



point in z since zw is different from zu, and 4 equivalent 
arriving points, 2 in z and 2 in u), 2 for zu and 1 for uw. 
The total symmetry factor is then 1, and the two inter- 
nal dashed lines are symmetrized together introducing a 
single normalization factor of ^. The graph in (|58|) has a 
symmetry factor of 1: there are 3-2 possibilities for xiy, 3 
for X2Z, 2-2 for the solid yz (the two lines in the loop are 
different and thus one must count also the starting legs), 
while the dashed yz is fixed. These exactly compensate 
the (1/2) (1/3!)^. In addition, there is one normalization 
factor of i. The same calculation holds for the graph 
([SO]) , while the one in (|56p has a symmetry factor of i 
(the two lines in the loop are indistinguishable) and no 
symmetrization. Some more examples are given in Fig. [1] 
In general, when the symmetrization of the internal 
dashed lines is normalized as in (|67p the total symmetry 
factor S of each A-th order graph is given by 



where P is the number of equivalent ways to place the 
lines. This is the same procedure to compute the sym- 
metry factors in standard QFT. However, unlike in the 
standard case, when counting the equivalent configura- 
tions one must distinguish solid and dashed lines. The 
resulting symmetry factor is therefore generically larger 
than the usual one, though still smaller than unity [and 
does not include all the factors of l/{n — 1)1 from the 
normalization of symmetrized dashed lines] . 

VII. SUMMARY AND CONCLUSIONS 

We are now able to summarize all the results obtained 
in the previous sections. An iterative solution of the 
equation of motion for the interacting operator field is 
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given, where higher order corrections are expressed in 
terms of retarded Green functions and lower order ones. 
The expectation value of n iterative solutions was shown 
to be exactly equivalent to an n-point in-in correlation 
function. On the other hand, the iterative solution pro- 
cedure suggested a Feynman-like picture of the various 
contributions, with a modified interpretation of the lines 
of the diagrams. 

Since two different kinds of line - solid or dashed - 
can be used, incquivalent versions of the same diagram 
can be generated simply by rearranging the disposition 
of the lines. Therefore, as a first step one has to draw all 
such possible versions, finding out all the allowed config- 
urations of solid and dashed lines. Secondly, the correct 
mathematical expression must be associated to each ele- 
ment of the graph. 

The prescriptions to build an allowed configuration are 
the following: 

• all solid lines form connected trees starting from an 
external point (there cannot be isolated solid lines) 

• all vertices are reached by at least one solid line 

• the total number of solid lines equals the number 
of vertices 

• there are no solid loops, nor connected trees reach- 
ing different external points 

• all the remaining lines are dashed 

The sum of all the allowed permutations of solid and 
dashed lines preserving the shape of the graph represents 
the total contribution of that Fcynman graph in the in-in 
formalism (see FiglT]). 

In turn, dashed lines belong to two different types: in- 
ternal (between vertices in the same tree) or external 
(connecting different trees). By definition, for any inter- 
nal dashed line there exists a unique loop having that 
as the only dashed line, which we have called the sym- 
mctrization loop. The lowest vertex of the loop, i.e. the 
one being the closest to the root of the tree, is the sym- 
mctrization vertex. In order to reconstruct the mathe- 
matical content of each graph, the following rules apply: 

• each vertex brings an integration factor —A j dya^ 

• each solid line corresponds to a retarded Green 
function Guiu, z) 

• an external dashed line represents the expectation 
value {(j)o{y)(f>o{z)) or {(j)o{z)(j)o{y)) . The order of 
the fields is fixed by the order of the respective 
trees in the correlation function 

• an internal dashed line not sharing the symmetriza- 
tion vertex contains the symmetrized expectation 

value i [{My)Mz)) + {Mz)My))] 

• several internal dashed lines sharing the same sym- 
metrization vertex arc symmetrized simultaneously 
using Eq. ^ or ^ 



• all graphs are multiplied by the symmetry factor 
(j68p . computed distinguishing between solid and 
dashed lines 

The above prescriptions introduce a useful practical 
method to reconnect the reformulation of the in-in for- 
malism given in Q with the Fcynman diagrams, pro- 
viding a quick and efficient tool to reconstruct all the 
perturbative terms without any calculation. The modi- 
fied Feynman rules describe loop correlation functions in 
terms of retarded Green functions and two-field expecta- 
tion values. This immediately refiects the decomposition 
outlined in Q of loop contributions into free fields and 
commutators of free fields. It thus allows to study in an 
easy way the late time behavior of these contributions, 
and analyze their dependence on powers of loga(t). A 
precise calculation is needed to investigate the possibil- 
ity of a resummation of all the powers of loga(<), which 
might eventually produce a dependence on some finite 
power of the scale factor. 

This modified diagrammatic formalism can also be 
used to study the divergences appearing in the loops. Al- 
though the integrals over time are basically understood, 
most diagrams are in fact plagued by divergences in the 
integrals over the three- momentum Q. A correct un- 
derstanding of these divergences is therefore necessary to 
investigate the behavior of the loop contributions, and 
to estimate their relevance with respect to the tree-level 
terms. The ultimate goal to achieve would be the removal 
of the divergences and the determination of a renormal- 
ized (at least up to some finite number of loops) effective 
potential. 

Finally, the analysis done here for a simple scalar field 
in a de Sitter background should be extended to a realis- 
tic case of the fluctuations of the inflaton field coupled to 
the metric. This case is of course more complicated be- 
cause of the presence of a higher number of fields and of 
the issue of gauge invariance, but the conceptual basis of 
the procedure outlined here and the setup of the new di- 
agrammatic formalism would not change. The resulting 
"renormalized" cosmological perturbation theory would 
be of greatest interest, since for practical purposes it 
might lead to different results. For example, the amount 
of non-Gaussianity hidden in loop contributions to the 
statistics of cosmological correlation functions might be 
larger than the standard prediction obtained with tree- 
level terms only. These aspects certainly deserve further 
investigation. 
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APPENDIX A: MORE ON ANTICOMMUTATORS 



As we have seen in Section llVi anticonimutators are 
introduced at each asymmetric vertex of an open graph, 
when the different branches have a different free field 
structure. Symmetric vertices do not need any anticom- 
mutator: either the free fields are all evaluated at the 
same time (as in final vertices) or they can be exchanged 
simply by renaming the integration variables. However, 
when open graphs are combined into closed graphs, the 
symmetry properties of each vertex also depend on what 
dashed lines are attached to it. 

As an example, let us consider the third order contribu- 
tion to the three-point function given by the contraction 
of the first term in with two free fields 4>o{x2) and 
4>oix3). This has the expectation value 



(Al) 



with no anticommutator in it (because we are free to 
exchange z and w). However, the result is different de- 
pending on which external free field is paired with which 
vertex. Actually, one gets 



4 (0o(^)0o(w)) {(f>Q{z)(l)o{x2)) (0o(w)0o(a;3)) 



(A2) 



-I- {(j)o{z)4)oix3)) {4'o{w)<j)o{x2)) ; 
after exchanging the variables, this becomes 

4{MZ)MX2)) {(l>oH4>3) {{^z^w) + {Mz)). (A3) 



where the internal dashed line is now symmetrized. How- 
ever, the symmetrization does not follow from the pres- 
ence of an anticommutator, but from the sum of the two 
different pairings with 00 (•'1^2) and (jj^ix^). 

The resulting closed graph representing this contribu- 
tion is 



Xl 



* X2 



(A4) 



even though in this case there is an intermediate vertex 
with no anticommutator associated, we can still follow 
the prescriptions outlined in Section |V| and interpret the 
internal dashed line as a symmetrized expectation value. 
Therefore, we do not need a different rule to symmetrize 
internal dashed lines within a completely symmetric sub- 
tree, even though the origin of the symmetrization is dif- 
ferent. 

We need anyway to pay attention to the symmetry fac- 
tor: if we consider the internal dashed line to be symmet- 
ric, there are 4 equivalent possibilities to pair the first ex- 
ternal leg and 2 for the second. This would overestimate 
the correct number of graphs given in (jA3p by a factor 
of 2. In general, we must add a factor ^ whenever the 
two solid branches to which the internal dashed line(s) 
is (are) attached are completely symmetric. This allows 
to adopt the same symmetrization rule for all internal 
dashed lines, and at the same time to use the standard 
combinatorial techniques when counting the number of 
equivalent graphs. 
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